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Motivated by the vast string landscape, we consider the shear viscosity to entropy density ratio 
in conformal field theories dual to Einstein gravity with curvature square corrections. After field 
redefinitions these theories reduce to Gauss-Bonnet gravity, which has special properties that allow 
us to compute the shear viscosity nonperturbatively in the Gauss-Bonnet coupling. By tuning of 
the coupling, the value of the shear viscosity to entropy density ratio can be adjusted to any positive 
value from infinity down to zero, thus violating the conjectured viscosity bound. At linear order in 
the coupling, we also check consistency of four different methods to calculate the shear viscosity, 
and we find that all of them agree. We search for possible pathologies associated with this class of 
theories violating the viscosity bound. 
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I. INTRODUCTION 



The AdS /conformal field theory (CFT) correspondence [J [3j S 13] has yielded many important insights into the 
dynamics of strongly coupled gauge theories. Among numerous results obtained so far, one of the most striking is the 
universality of the ratio of the shear viscosity t] to the entropy density s d, H, 0, [1] 

- = ^ (1-1) 

s An ^ ' 

for all gauge theories with an Einstein gravity dual in the limit N oo and A — s- oo. Here, N is the number of 
colors and A is the 't Hooft coupling. It was further conjectured in |8l| that (jl.ip is a universal lower bound [the 
Kovtun-Starinets-Son (KSS) bound] for all materials. So far, all known substances including water and liquid helium 
satisfy the bound. The systems coming closest to the bound include the quark-gluon plasma created at Relativistic 
Heavy Ion Colhder (RHIC) [§, [13) [HI 112, [Hj [13] and certain cold atomic gases in the unitarity limit (see e.g. [11]). v/s 
for pure gluon QCD slightly above the deconfinement t emp erature has also been calculated on the lattice recently [ig] 



and is about 30% larger than dTT]). See also pj]. See [18|, [ij, l20|, |2l|, 122] for other discussions of the bound. 

Now, as stated above, the ratio (jl.ip was obtained for a class of gauge theories whose holographic duals are dictated 
by classical Einstein gravity (coupled to matter). More generally, string theory (or any quantum theory of gravity) 
contains higher derivative corrections from stringy or quantum effects, inclusion of which will modify the ratio. In 
terms of gauge theories, such modifications correspond to 1/A or 1/N corrections. As a concrete example, let us 
take Af — 4: super- Yang-Mills theory, whose dual corresponds to type IIB string theory on AdS^ x S^. The leading 
order correction in 1/ A arises from stringy corrections to the low-energy effective action of type IIB supergravity, 
schematically of the form a'^R^. The correction to 77/s due to such a term was calculated in [H, [l^]. It was found 
that the correction is positive, consistent with the conjectured bound. 

In this paper, instead of limiting ourselves to specific known string theory corrections, we explore the modification of 
77/5 due to generic higher derivative terms in the holographic gravity dual. The reason is partly pragmatic: other than 
in a few maximally supersymmetric circumstances, very little is known about forms of higher derivative corrections 
generated in string theory. Given the vastness of the string landscape [2^ , one expects that generic corrections do 
occur. Restricting to the gravity sector in AdS^, the leading order higher derivative corrections can be written as^ 

/ = j d^x ^ {R - 2A + L2 (aii?2 ^ a^R^.R^^" + a^R'^'^P'' R,,,p,)) , (1.2) 

where A = — -p- and for now we assume that ^ -p- ^ 1. Other terms with additional derivatives or factors of R 
are naturally suppressed by higher powers of j^. String loop (quantum) corrections can also generate such terms, but 
they are suppressed by powers of gs and we will consistently neglect them by taking —> limit. ^ To lowest order 
in ai the correction to rj/s will be a linear combination of a^'s, and the viscosity bound is then violated for one side 
of the half-plane. Specifically, we will find 

- = ^(l"8a3) + 0(a?) (1.3) 
s 47r 

and hence the bound is violated for > 0. Note that the above expression is independent of ai and a2. This can 
be inferred from a field redefinition argument (see Sec lIIC|) . 
How do we interpret these violations? Possible scenarios are: 

1. The bound can be violated. For example, this scenario would be realized if one explicitly finds a well-defined 
string theory on AdS^ which generates a stringy correction with > 0. (See (27| for a plausible counterexample 
to the KSS bound.) 

2. The bound is correct (for example, if one can prove it using a field theoretical method), and a bulk gravity 
theory with as > cannot have a well-defined boundary CFT dual. 



^ Our conventions are those of [26ll . In this section we suppress Gibbons-Hawking surface terms. 

^ Note that to calculate corrections, all the light fields must be taken into account. In addition, the calculation of rj/s could be more 
subtle once we begin to include quantum effects. 
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(a) The bulk theory is manifestly inconsistent as an effective theory. For example, it could violate bulk causality 
or unitarity. 

(b) It is impossible to generate such a low-energy effective classical action from a consistent quantum theory 
of gravity. In modern language we say that the theory lies in the swampland of string theory. 

Any of these alternatives, if realized, is interesting. Needless to say, possibility 1 would be interesting. Given that 
recent analyses from RHIC data Elj [Hj [Hj [3 indicate the rj/s is close to (and could be even smaller than) the 
bound, this further motivates to investigate the universality of the KSS bound in holographic models. 

Possibility 2(a) should help clarify the physical origin of the bound by correlating bulk pathologies and the violation 
of the bound. Possibility 2(b) could provide powerful tools for constraining possible higher derivative corrections in 
the string landscape. Note that while there are some nice no-go theorems which rule out classes of nongravitational 
effective field theories [11] (also see [2§|), the generalization of the arguments of [111 to gravitational theories is subtle 
and difficult. Thus, constraints from AdS/CFT based on the consistency of the boundary theory would be valuable. 

In investigating the scenarios above, Gauss-Bonnet (GB) gravity will provide a useful model. Gauss-Bonnet gravity, 
defined by the classical action of the form [30[ 



/ = 



1 



i?-2A+-'^^^^2,o2 



(1.4) 



has many nice properties that are absent for theories with more general ratios of the a^'s. For example, expanding 
around flat Minkowski space, the metric fluctuations have exactly the same quadratic kinetic terms as those in Einstein 
gravity. All higher derivative terms cancel ^Q] . Similarly, expanding around the AdS black brane geometry, which will 
be the main focus of the paper, there are also only second derivatives on the metric fluctuations. Thus small metric 
fluctuations can be quantized for finite values of the parameter Xgb-^ Furthermore, crucial for our investigation is 
its remarkable feature of solvability: sets of exact solutions to the classica l eq uation of motion have been obtained 
[31I . [3^ and the exact form of the Gibbons-Hawking surface term is known [33|. 

Given these nice features of Gauss-Bonnet gravity, we will venture outside the regime of the perturbatively corrected 
Einstein gravity and study the theory with finite values of Xqb- To physically motivate this, one could envision that 
somewhere in the string landscape Xgb is large but all the other higher derivative corrections are small. One of the 
main results of the paper is a value of ri/s for the CFT dual of Gauss-Bonnet gravity, nonperturbative in Xgb-'^ 



^ = ]-[1-4Xgb] 
s 47r 



(1.5) 



We emphasize that this is not just a linearly corrected value. In particular, the viscosity bound is badly violated as 
Xgb — > J- As we will discuss shortly, Xgb is bounded above by j for the theory to have a boundary CFT, and r//s 
never decreases beyond 0. 

Given the result (jl.Sp for Gauss-Bonnet, if the possibility 2(a) were correct, we would expect that pathologies would 
become easier to discern in the limit where 77/s — > 0. We will investigate this line of thought in Sec lIVI On the other 
hand, thinking along the line of possibility 1, the Gauss-Bonnet theory with Xgb arbitrarily close to j may have a 
concrete realization in the string landscape. In this case, there exists no lower bound for rj/s, and investigating the 
CFT dual of Gauss-Bonnet theory should clarify how to evade the heuristic mean free path argument for the existence 
of the lower bound (presented in, e.g., Q). 

The plan of the paper is as follows. In SecHH we review various properties of two-point correlation functions and 
outline the real-time AdS/CFT calculation of the shear viscosity. We then explicitly calculate the shear viscosity 
for Gauss-Bonnet theory in Sec lIIII In Sec lIVl we seek possible pathologies associated with theories violating the 
viscosity bound. There, we will find a curious new metastable state for large enough Xgb- Finally in Sec|Vl we 
conclude with various remarks and speculations. To make the paper fairly self-contained, various appendices are 
added. In particular, quasinormal mode calculations of the shear viscosity are presented in Appendix [B] and one using 
the membrane paradigm in Appendix [D] 



^ Generic theories in 111.21 1 contain four derivatives and a consistent quantization is not possible other than treating higher derivative terms 
as perturbations. 

We have also computed the value of r)/s for Gauss-Bonnet gravity for any spacetime dimension D and the expression is given in US. 261 1 . 
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II. SHEAR VISCOSITY IN THEORIES: PRELIMINARIES 
A. Two-point correlation functions and viscosity 

Let us begin by collecting various properties of two-point correlation functions, following [s^jIUjllBl (see also [37|). 
Consider retarded two-point correlation functions of the stress energy tensor T^^, of a CFT in 3 -|- 1-dimensional 
Minkowski space at a finite temperature T 

G^.M^,q) = -t [ dtdxe''^'-'^^-'^e{t){[T^,{t,x),T^0{O,O)]). (2.1) 



They describe linear responses of the system to small disturbances. It turns out that various components of (|2.ip 
can be expressed in terms of three independent scalar functions. For example, if we take spatial momentum to be 
q (0, 0, q), then 

1 1 cj^ 2 Lu'^ 
Gi2,i2 = -G^{uj,q), Gi3,i3 = -^ -Gi{uj,q), 633,33 = --^ —G2{uJ,q), (2.2) 

and so on. Kt q — Q all three function Gi_2.3('^, 0) are equal to one another as a consequence of rotational symmetry. 

When LL!,\q\ ^ T one expects the CFT plasma to be described by hydrodynamics. The scalar functions Gi,2,3 
encode the hydrodynamic behavior of shear, sound, and transverse modes, respectively. More explicitly, they have 
the following properties: 

• Gi has a simple diffusion pole dX to ~ —iDq^, where 

D = ^ = ^1 (2.3) 

e + P Ts ^ ' 

with e and s being the energy and entropy density, and P the pressure of the gauge theory plasma. 

• G2 has a simple pole ai lu — ±Csq~ iT sq^ , where Cg is the speed of sound and F^. is the sound damping constant, 
given by (for conformal theories) 

' 3Ts ^ ' ' 

• rj can also be obtained from Gi.2,3 at zero spatial momentum by the Kubo formula, e.g., 

77== lim -ImGi2i2(t^,0) (2.5) 

Equations (|2.3p - (|2.5p provide three independent ways of extracting 'q/s. We provide calculations utilizing the first 
two in Appendix [B] A calculation utilizing the Kubo formula (12. 5p is easier, and we will explicitly implement it for 
Gauss-Bonnet theory in Sec lIIII In the next subsection, we outline how to obtain retarded two-point functions within 
the framework of the real-time AdS / CFT correspondence. 



B. AdS/CFT calculation of shear viscosity: Outline 

The stress tensor correlators for a boundary CFT described by p.2p or (|1.4p . can be computed from gravity as 
follows. One first finds a black brane solution (i.e. a black hole whose horizon is R'^) to the equations of motion of 
(II. 2p or (|1.4p . Such a solution describes the boundary theory on R^'^ at a temperature T, which can be identified with 
the Hawking temperature of the black brane. The entropy and energy density of the boundary theory are given by 
the corresponding quantities of the black brane. The fluctuations of the boundary theory stress tensor are described 
in the gravity language by small metric fluctuations /i^jy around the black brane solution. In particular, after taking 
into account of various symmetries and gauge degrees of freedom, the metric fluctuations can be combined into three 
independent scalar fields 4>a,a = 1,2, 3, which are dual to the three functions Ga of the boundary theory. 

To find Go, one could first work out the bulk two-point retarded function for (j)a and then take both points to the 
boundary of the black brane geometry. In practice it is often more convenient to use the prescription proposed in [38| , 
which can be derived from the real-time AdS/CFT correspondence [s^. Let us briefly review it here: 



5 



1. Solve the linearized equation of motion for 4>aif", k) with the following boundary conditions: 

(a) Impose the infalling boundary condition at the horizon. In other words, modes with timelike momenta 
should be falling into the horizon and modes with spacelike momenta should be regular. 

(b) Take r to be the radial direction of the black brane geometry with the boundary at r = oo. Require 

Mr;k)l^i^Ja{k), k^{LO,q), (2.6) 

where e ^ imposes an infrared cutoff near the infinity of the spacetime and Ja{k) is an infinitesimal 
boundary source for the bulk field ipaii", k). 

2. Plug in the above solution into the action, expanded to quadratic order in 4)a{'r\ k). It will reduce to pure surface 
contributions. The prescription instructs us to pick up only the contribution from the boundary aX r = -. The 
resulting action can be written as 

S^~IJ M-k)Mk,r)Jaik)l^^ . (2.7) 

Finally the retarded function Ga{k) in momentum space for the boundary field dual to 0a is given by 

Ga{k) = \imTa{,k,r) ^ . (2.8) 



Using the Kubo formula (j2.5p . we can get the shear viscosity by studying a mode (j)^ with q — Q in the low- frequency 
limit Lo Q. We will do so in the next section. Alternatively, using (|2.3|) or \2A\ . we can read off the viscosity from 
pole structures of retarded two-point functions. Such a calculation is a bit more involved and will be performed in 
Appendix [B] 

The above prescription for computing retarded functions in AdS/CFT works well if the bulk scalar field has only two 
derivatives as in Gauss-Bonnet case ()1.4p . If the bulk action contains more than two derivatives, complications could 
arise even if one treats the higher derivative parts as perturbations. For example, one needs to add Gibbons-Hawking 
surface terms to ensure a well-defined variational problem. A systematic prescription for doing so is, however, not 
available at the moment beyond the linear order. Thus there are potential ambiguities in implementing (j2.8p .^ Clearly 
these are important questions which should be explored more systematically. At the level, as we describe below 
in Sec lII C[ all of our calculations can be reduced to the Gauss-Bonnet case in which these potential complications do 
not arise. 



C. Field redefinitions in tlieories 



We now show that to linear order in a^, rj/ s for (jl.2p is independent of ai and a2- It is well known that to linear 
order in ai, one can make a field redefinition to remove the and Rf^^R^'^ term in p.2p . More explicitly, in (|1.2p 
set as = and take 

g^^i' = fifj^i' + a2L^R^_iu - -^("2 + 2ai)gf_,^R, (2.9) 
where R denotes the Ricci scalar for g^^i, and so on. Then (jl.2p becomes 

^ ^ TTtV / + - 2A) + 0{a') = / V^{R ~ 2A) + 0{a^) (2.10) 

with 

^ 2AL2 ^ ^ ~ A , , 

/C = ^ (5ai + a2) , A=YT^- (2.11) 



^ In [24|. such additional terms do not appear to affect tlie calculation at the order under discussion there. 
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It follows from (12. 9p that a background solution g^"^ to (|1.2p (with = 0) is related to a solution g^"^ to (|2.10p by 

- 2 

ds^^A^dso, A = l-— . (2.12) 

o 

The scaling in (|2.12p does not change the background Hawking temperature. The diffusion pole (|2.3p calculated using 
(|2.10p around then gives the standard result D = [34]. Thus we conclude that rj/s = for (|1.2p with 
ck3 = 0. Then to linear order in a^, 77/s can only depend on a^. To find this dependence, it is convenient to work 
with the Gauss-Bonnet theory (II. 4p . Gauss-Bonnet gravity is not only much simpler than ()1.2p with generic (23 ^ 0, 
but also contains only second derivative terms in the equations of motion for hp^i, , making the extraction of boundary 
correlators unambiguous. 



III. SHEAR VISCOSITY FOR GAUSS-BONNET GRAVITY 



In this section, after briefly reviewing the thermodynamic properties of the black brane solution, we compute the 
shear viscosity for Gauss-Bonnet gravity (|1.4p nonperturbatively in Xgb- Here, we follow the outline presented in the 
previous section, with the Kubo formula p.Sp in mind. In AppendixlBj we extract rj/s from the shear channel (|2.3p and 
the sound channel (|2.4p (perturbatively in Xgb)- There we also find that the sound velocity remains at the conformal 
value = ^ as it should. In Appendix [D] we provide a membrane paradigm calculation, again nonperturbatively in 
Xgb- All four methods give the same result. 



A. Black brane geometry and thermodynamics 



Exact solutions and thermodynamic properties of black objects in Gauss-Bonnet gravity (|1.4p were discussed in |32l | 
(see also [13, ST], SI])- Here we summarize some features relevant for our discussion below. The black brane 
solution can be written as 




where 



fir) 



1 



^Xgb 



4A 



GB 



In p.ip , is an arbitrary constant which specifies the speed of light of the boundary theory. Note that as r 

-j^, with = -[l + ^/l- AXgb ] 



fir) 



(3.1) 



(3.2) 



(3.3) 



It is straightforward to see that the AdS curvature scale of these geometries is aL.^ If we choose — a, then the 
boundary speed of light is unity. However, we will leave it unspecified in the following. We assume that Xgb < j- 
Beyond this point, p.4p does not admit a vacuum AdS solution, and cannot have a boundary CFT dual. In passing, 

we note that while the curvature singularity occurs at r = for Xgb > 0, it shifts to r = r+ ^1 — 4^^^ for 
Xgb < 0. 

The horizon is located at r = r_|_ and the Hawking temperature, entropy density, and energy density of the black 
brane are ^ 



r+ 



(3.4) 



^ Here we note that the Gauss-Bonnet theory also admits another background with the curvature scale a L where a? = ^ (l — \/l — i^Gs) ■ 
Even though this remains an asymptotically AdS solution for Xqb > 0, we do not consider it here because this background is unstable 
and contains ghosts |31j . 

Note that for planar black branes in Gauss-Bonnet theory, the area law for entropy still holds 44] . This is not the case for more general 
higher-derivative-corrected black objects. 
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If we fix the boundary theory temperature T and the speed of fight to be unity (taking iVj = a), the entropy and 
energy density are monotonicaUy increasing functions of Agb, reaching a maximum at Xgb — \ and going to zero as 
Acs -oo. 

To make our discussion self-contained, in Appendix |^ we compute the free energy of the black brane and derive 
the entropy density. In particular, we show that the contribution from the Gibbons-Hawking surface term to the free 
energy vanishes. 



4G 



B. Action and equation of motion for the scalar channel 

To compute the shear viscosity, we now study small metric fiuctuations (j) — h^2 around tfie black brane background 
of the form 

ds^ = -f{r)N^dt^ + Jf^dr^ + J2 {TI'^'^^^ + 2(j){t,x,r)dxidx^ . (3.6) 

We will take (j) to be independent of xi and X2 and write 

dkjdq 

W) 

For notational convenience, let us introduce 



<f>{t,x,r)^ I j^cj,ir;k)e~''^'+^'^^\ fc = (c., 0, 0, g), (/.(r; -fc) = ^ fc) • (3.7) 



z^-, q^-q. /=;:r/-2A^|l-\/l-4AGB + -^). (3.8) 



Then, at quadratic order, the action for ^ can be written as 

^ = / ^^^^^^ 6'(fci, fc2) with 
J (27r)^ 

5(fci = 0, fc2 = 0) = 1 dz'^ {Kid.ct^f - K24>^ + d.{K,^^)) , (3.9) 



where 



and (fy^ should be understood as a shorthand notation for 0(z; k)(j){z, —k). Here, S is the sum of the bulk action (|1.4[) 
and the associated Gibbons-Hawking surface term [s^ . The explicit expression for K3 will not be important for our 
subsequent discussion. 

The equation of motion following from (|3.9[) is^ 

K(j>" + K'(I)' + K2(t>^Q , (3.11) 

where primes indicate partial derivatives with respect to z. Using the equation of motion, the action (j3.9p reduces to 
the surface contributions as advertised in Sec lIIBl 



Sik, = 0, fc2 = 0) = -ic J ^ {K^'cb + K^ict^^) Surface . (3.12) 



The prescription described in Sec lII Bl instructs us to pick up the contribution from the boundary at z ^ 00. Here, 
the term proportional to K'i will give rise to a real divergent contact term, which is discarded. 



An easy way to get the quadratic action 113.91 1 is to first obtain the linearized equation of motion and then read off K and if 2 from it. 
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A curious thing about (|3.9p is that for all values of z, both K and K2 (but not K^i) are proportional to \ — Xgb-^ 
Thus other than the boundary term the whole action p.9p vanishes identically at Acs = Nevertheless, the equation 
of motion (|3.1ip remains nontrivial in the limit Xgb ^ 3 as the j — Xgb factor cancels out. Note that the correlation 
function does not necessarily go to zero in this limit since it also depends on the behavior of the solution to (|3.1ip 
and the limiting procedure (|3.12p . As we will see momentarily, as least in the small frequency limit it does become 
zero with a vanishing shear viscosity. 



C. Low- frequency expansion and the viscosity 

General solutions to the equation of motion (|3.1ip can be written as 

k) = ain{k)<i)in(z\ k) + aout{k)4>out{z] k) , (3.14) 

where and (j)out satisfy infalling and outgoing boundary conditions at the horizon, respectively. They are complex 
conjugates of each other, and we normalize them by requiring them to approach 1 as z ^ 00. Then, the prescription 
of Sec lII Bl corresponds to setting 

a-inik) = J{k) , aout{k) ^ , (3.15) 

where J(fc) is an infinitesimal boundary source for the bulk field (f>. 

More explicitly, as z — * 1, various functions in p. lip have the following behavior 

It follows that near the horizon z = 1, equation ()3.1ip can be solved by (for q = 0) 

<j}{z) ^ {z - if ^ {z - 1)^T^ (3.17) 

with the infalling boundary condition corresponding to the negative sign. To solve (|3.1ip in the small frequency limit, 
it is convenient to write 

0.„(z; k) = e"^(^)'"(^) (^1 - «^5i(^) + 0{Cj^f)^ , (3.18) 

where we require gi{z) to be nonsingular at the horizon z = 1. We show in Appendix [Cl that gi is a nonsingular 
function with the large z expansion 



Therefore, with our boundary conditions (|3.15p . we find 



{z;k) = J{k) 



(3.20) 



This is the right asymptotic behavior for the bulk field (j) describing metric fluctuations since the CFT stress tensor 
has conformal dimension 4. 

Plugging (|3.20p into (|3.12p and using the expressions for C and K in (|3.10p . the prescription described in Sec lII Bl 
gives 

ImGi2.i2(t^, 0) = oj—^ f ^) (1 - 4Ags) + 0{uj^). (3.21) 
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This can be seen by using the following equation in K and K2 



z^-2XGBf 
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Then, the Kubo formula (|2.5p yields 



Finally, taking the ratio of (|3.22p and p.Sp we find that 

1 = ^{1-4Xgb). (3.23) 

This is nonperturbative in Acs- Especially, the linear correction is the only nonvanishing term.^" 
We now conclude this section with various remarks: 

1. Based on the field redefinition argument presented in Sec lII CI one finds from (|3.23p that for (|1.2|) . 

- = ^ (1 - Sas) + 0(a?). (3.24) 
s 47r 

We have also performed an independent calculation of rj/s (without using field redefinitions) for (|1.2p using all 
three methods outlined in Sec lII Al and confirmed (|3.24p . 



2. The ratio rj/s dips below the viscosity bound for Acs > in Gauss-Bonnet gravity and for as > in (|1.2p . 
In particular, the shear viscosity approaches zero as Xqe ^ j for Gauss-Bonnet. Note that the whole off-shell 
action becomes zero in this limit. It is likely the on-shell action also vanishes, implying that the correlation 
function could become identically zero in this limit. 

3. Fixing the temperature T and the boundary speed of light to be unity, as we take Xqb — > — oo, 77 ~ {—Xgb)^ — > 
00. In contrast the entropy density decreases as s '--^ (—Xgb)^^ ~* 0. 

4. The shear viscosity of the boundary conformal field theory is associated with absorption of transverse modes by 
the black brane in the bulk. This is a natural picture since the shear viscosity measures the dissipation rate of 
those fiuctuations: the quicker the black brane absorbs them, the higher the dissipation rate will be. For example, 
as Xgb ^ —00, rj/s approaches infinity; this describes a situation where every bit of the black brane horizon 



devours the transverse fluctuations very quickly. In this limit the curvature singularity at z = (^1 — ^^ge 

approaches the horizon and the tidal force near the horizon becomes strong. On the other hand, as Xgb ^ j, 
Tj/s —> and the black brane very slowly absorbs transverse modes. 

5. The calculation leading to (|3.23p can be generalized to general D spacetime dimensions and one finds for 
L> > 4 + 



rj 1 

s An 



(D-l), 
(15-3) 



(3.26) 



Here again Xgb is bounded above by j. Thus for £> > 4 -I- 1, 77 never approaches zero within Gauss-Bonnet 
theory. For _D = 3 + lor2-|-l, in which case the Gauss-Bonnet term is topological, there is no correction to 
r]/s. 

6. In Appendix [Dl we obtain the same result (|3.23p using the membrane paradigm [6i]. Thus when embedded 
into the AdS/CFT correspondence, the membrane paradigm correctly captures the infrared (hydrodynamic) 
sector of the boundary thermal field theory. Further, we see something interesting in its derivation. There, the 
diffusion constant is expressed as the product of a factor evaluated at the horizon (jD22p and an integral from 
the horizon to infinity (|D23p . In the limit Xgb ~> -j, it is the former that approaches zero. 



It would be interesting to find an explanation for vanishing of higher order corrections. 

We note that for Xqb = ^ in 4 -(- 1 spacetime dimension, the radial direction of the background geometry resembles a Bafiados- 

Teitelboim-Zanelli (BTZ) black brane. 

For general dimensions we use the convention 

S = —^—r- f d°x [R-2K + aGBL^iR" - iRf.uR'"' + iJ^.pai?"""")] (3.25) 
levrGjv J 

with A = -L£^:M^^ and Xgb = {D - 3){D - A)o,gb- 
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FIG. 1: c^iz) (vertical axis) as a function of z (horizontal axis) for Xgb = 0.08 (left panel) and Xqb =0.1 (right panel). For 
Xgb < yIq , Cg is a monotonically increasing function of z. When Xgb > y|q, as one decreases z from infinity, Cj increases from 
1 to a maximum value at some z > 1 and then decreases to as 2; —> 1 (horizon). 



IV. CAUSALITY IN BULK AND ON BOUNDARY 



In this section we investigate if there are causaUty problems in the bound-violating theories discussed above. First 
we will discuss the bulk causal structure. Then we discuss a curious high-momentum metastable state in the bulk 
graviton wave equation that may have consequences for boundary causality. The analysis in this section is refined 
in psj where we indeed see a precise signal of causality violation for Ag b > y§q ■ 



A. Graviton cone tipping 



As a consequence of higher derivative terms in the gravity action, graviton wave packets in general do not prop- 
agate on the light cone of a given background geometry. For example, when Xqe 7^ 0, the equation (|3.1ip for the 
propagation of a transverse graviton differs from that of a minimally coupled massless scalar field propagating in the 
same background geometry (|3.ip . To make the discussion precise, let us write (we will consider only xi^2-independent 
waves) 

0(i,r,X3) = e-*"*+*-'^+''?^^(/.e„(t,r,X3). (4.1) 



Here, (j)en is a slowly-varying envelope function, and we take the limit k — (oj, fc^, 0, 0, g) — > 00. In this limit, the 
equation of motion (|3.1ip reduces to 

k'^k'-gf, « 0, (4.2) 

where 

dsls = gfjx^dx^ = f{r)Nl (^-d^ + ^dx'^ + -^dr\ (4.3) 



In (OD 



\ _ ^GBJ 2 



z 

can be interpreted as the local "speed of graviton" on a constant r-hypersurface. = ^'f'j^'' introduced in the 
second equality in (|4.4[) is the local speed of light as defined by the background metric (|3.ip . Thus the graviton cone 
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in general does not coincide with the standard null cone or light cone defined by the background metric. -'^^ A few 
more comments about graviton cone are found at the end of Appendix iDl 

In the nongravitational boundary theory there is an invariant notion of light cone and causality. At a heuristic 
level, a graviton wave packet moving at speed Cg{z) in the bulk should translate into disturbances of the stress tensor 
propagating with the same velocity in the boundary theory. It is thus instructive to compare Cg and Cb with the 
boundary speed of light, which we now set to unity by taking — a {a was defined in ()3.3|1 ). At the boundary 
(z = oo) one finds that Cg{z) = Cb{z) = 1. In the bulk, the background local speed of light Cb is always smaller than 1, 
which is related to the redshift of the black hole geometry. The local speed of graviton Cg(z), however, can be greater 
than 1 for certain range of z if Xgb is sufficiently large. To see this, we can examine the behavior of near z = oo, 

Cg{z) ~ ^ ^ -i + 0{z ), z^oo, biiXoB)^ 2(1 -4Agb) ' 

biiXas) becomes positive and thus increases above 1 if Xqb > y§o' ^'^'^ such a Xgb, as we decrease z from infinity, 
Cg will increase from 1 to a maximum at some value of z and then decrease to zero at the horizon. See Fig. [T] for 
the plot of Cg(z) as a function of z for two values of Xgb- When Xgb = y§o finds that the next order term in 
(|4.6|) is negative and thus Cg does not go above 1. Also note that Xgb \, &i(Agb) goes to plus infinity. Thus 
heuristically, in the boundary theory there is a potential for superluminal propagation of disturbances of the stress 
tensor. 

In (45I ] we explore whether such bulk graviton cone behavior can lead to boundary causality violation by studying the 
behavior of graviton null geodesies in the effective geometry. There, we indeed see causality violation for Xgb > j^H' 



B. New metastable states at high momenta {Xgb > y§q) 

We now study the behavior of the full graviton wave equation. Let us recast the equation (|3.1ip in Schrodinger 
form. For this purpose, we introduce 



Then (jS.lip becomes 



with 




(4.7) 



(4.8) 



Viy) ^ f4iz) + Vi, yi(y) = 1|- ^ \^B" + ^B'] , (4.9) 

where Cg(z) was defined in (|4.4p . The advantage of using (|4.8p is that qualitative features of the full graviton 
propagation (including the radial direction) can be inferred from the potential V{y), since we have intuition for 
solutions of the Schrodinger equation. Since y is a monotonic function of z, below we will use the two coordinates 
interchangeably in describing the qualitative behavior of V{y). 

One can check that Vi(z) is a monotonically increasing function for any Xgb > (note Vi{z) +00 as z 00). 
For Xgb < jin, c'i{z) is also a monotonically increasing function as we discussed in the last subsection and the whole 



Note that 



' _ 1 - XcBf' ^ 1 - 4Agb + 12^ 
■2 l^Acs/i 1--4Ags+4^ 



s 



(4.5) 

b ^ 

and in particular the ratio is greater than 1 for Xqe > 0- Note that bulk causality and the existence of a well-posed Cauchy problem 
do not crucially depend on reference metric light cones and such tipping is not a definitive sign of causality problems. Also for any 
value of Xgb, the graviton cone coincides with the light cone in the radial direction. If not, we could have argued for the violation of 
the second law of thermodynamics following l46l |47|| . Further note that for Xgb < ~§i there exists a region outside the horizon where 
Cg<0 which will lead to the appearance of tachyonic modes, following [48|. We have not explored the full significance of this instability 
here since it is not correlated with the viscosity bound. 

In fact coefficients of all higher order terms in I/2 expansion become divergent in this limit. 



12 




FIG. 2: V{z) - (vertical axis) as a function of z (horizontal axis) for Xgb = 0.08 and q = 500 (left panel) and for Xgb = 0.1 
and q = 500 (right panel). V{z) is a monotonically increasing function of z for Xgb < j§g, but develops a local minimum for 
Xgb > Y§Q with large enough q. 

V{z) is monotonic. When Xgb > y§q , there exists a range of z where Cg(z) decreases with increasing z for sufficiently 
large z. Thus V{z) can now have a local minimum for sufficiently large q. For illustration, see Fig. [2] for the plot of 
V(z) as a function z for two values of Xgb- 

Generically, a graviton wave packet will fall into the black brane very quickly, within the time scale of the inverse 
temperature ^ (since this is the only scale in the boundary theory). Here, however, precisely when the local speed of 
graviton Cg can exceed 1 (i.e. for Xgb > y^q)' ^(•^) develops a local minimum for large enough q and the Schrodinger 
equation (14. 8p can have metastable states living around the minimum. Their lifetime is determined by the tunneling 
rate through the barrier which separates the minimum from the horizon. For very large q this barrier becomes very 
high and an associated metastable state has lifetime parametrically larger than the timescale set by the temperature. 
In the boundary theory, these metastable states translate into poles of the retarded Green function for T^y in the 
lower half-plane. The imaginary part of such a pole is given by the tunneling rate of the corresponding metastable 
state. Thus for Xgb > y§qi in boundary theory we find new quasiparticles at high momenta with a small imaginary 
part.^^ 

In [4^ , we confirm that those long-lived quasiparticles give rise to causality violation for Xg b > • 

V. DISCUSSION 

In this paper we have computed rj/s for Gauss-Bonnet gravity using a variety of techniques. We have found that 
the viscosity bound is violated for Xgb > and have looked for pathologies correlated to this violation. For small 
positive Xgb we have not found any. The violation of the bound becomes extreme as Xgb — > j where 77 vanishes. 
We have focused our attention on this region to find what unusual properties of the boundary theory could yield a 
violation not only of the bound but also of the qualitative intuitions suggesting a lower bound on rj/s. Above we also 
have discussed a novel quasiparticle excitation. In [i^, causality violation is firmly established for Xgb > 

It is also instructive to examine the behavior of the zero temperature theory as Xgb ^ \- Basic parameters 
describing the boundary CFT are the coefficients of the 4D Euler and Weyl densities called a and c respectively. 
These have been computed first in [5^ , and for Gauss-Bonnet gravity in [5l| . Their results indicate that 

c - {1-4:Xgb)K 

a ^ (3(1-4Agb)^ -2). (5.1) 

The parameter c is related to the two-point function of a boundary stress tensor which is forced by unitarity to be 
positive. (|5.ip shows that c vanishes at Xgb = j demonstrating the sickness of this point. For Ags a bit less than j 
the stress tensor couples very weakly in a system with a large number of degrees of freedom. This is peculiar indeed. 
In the bulk it seems that gravity is becoming strongly coupled there. 



A similar type of long-lived quasiparticles exist for Af = A S YM theory on [4S| , but not on R,^ . 
This can also be seen from the derivations in Sec lIIII 
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The coefficient a vanishes at Xqb = The significance of this is unclear. 

More generally, we believe it would be valuable to explore how generic higher derivative corrections modify various 
gauge theory observables. This is important not only for seeing how reliable it is to use the infinite 't Hooft coupling 
approximation for questions relevant to QCD, but also for achieving a more balanced conceptual picture of the 
strong coupling dynamics. Furthermore, this may generate new effective tools for separating the swampland from the 
landscape. 

As a cautionary note we should mention that pathologies in the boundary theory in regions that violate the viscosity 
bound may not be visible in gravitational correlators, at least when gs = 0. As an example consider the a'^R^ terms 
discussed in [23|- For positive a' , the physical case, the viscosity bound is preserved. But the bulk effective action 
can equally be studied for a' negative. Here gravitational correlators can be computed and will violate the viscosity 
bound. The only indication of trouble in the boundary theory at = will come from correlators of string scale 
massive states, whose mass and CFT conformal weight ~ l/{a')2 ^ an imaginary number! 
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APPENDIX A: THERMODYNAMIC PROPERTIES OF GB BLACK HOLES 



1. Free energy 



It is easy to confirm that the following metric is a stationary point of the Gauss-Bonnet action (|1.4I 



i-dr^ + ^ f y dxU with 



/( 

/(.) = (l - yi-4Ac.+4A,.(!±)^j . (Al) 

First note that the Hawking temperature is 

r(r+) = = iVjl^. (A2) 

ZTT y^grr dr TT 

To get the free energy F[T] of the macroscopic configuration (|A1[) . we note the following correspondence in the 
classical limit: 

e-TnT]^Z[T]^e-'^^l (A3) 

Here, I[T] is the Euclidean action of the configuration with temperature T. Evaluating the Euclideanized bulk action 
for Gauss-Bonnet gravity (|1.4p with the background metric (|Aip . we find 

/buik[n^+)] = "I^^/ j d^x,^[R-2K+^L\R^ -AR^.^R'^'' + R^,p,W'P'')] 
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167rGAr T L5 Acs 



(i2Agb - 5 + 5^1 -4AGi3) - 4Ags 



2Agb 



^/l~^4A^ 



(A4) 



We regulate this result by subtracting the Euclidean action of the AcB-modified pure AdS space (obtained by setting 
r+ = in (Ml ) 



1 



ri 1 



16ttGn^^T' L5 Ags 



(i2Agb-5 + 5v/1-4Agb 



with T'{T) chosen so that the geometries at r = Tmax agree (J]. Quantitatively, 



(A5) 



r2 

max 



L2 2AGi3 



1 - Wl -4Agb +4AGi3- 



(A6) 



All in all, we get 



f|Ti = T(7t,.,i.ir]-C;jrmi) 



The entropy density is then given by 



(A7) 



(A8) 



2. Vanishing of Gibbons-Hawking contribution 



To be complete, we need to show that there is no contribution to the free energy from the Gibbons-Hawking surface 
term when we regulate with the background subtraction method presented above. This can be shown explicitly. For 
the black brane solution, the Gibbons-Hawking contribution is^^: 



lGH[T(r+)] = 



1 



-V, 



16ttGj\[ 

1 ^3^^(_.,3yTT«^ 



6Ags <3L^{^ 



drj 

r 



r^4 

max 



1 - yi - 4AGi3 

Agb 



4:TtGn 

A similar expression is obtained for pure AdS space. With the choice (jA6p . we obtain 

Igh[T] 



/grr(T)]^o 



1 



VI - 4Agb 



(A9) 



(AlO) 



APPENDIX B: 7?/s FROM SHEAR AND SOUND CHANNEL POLES 



Our calculation in this appendix follows the techniques developed in [36l |. 

Consider a perturbation of the background metric of the form h^^, — ft,^i,(r)e~'"*"'"*'^^, with ^,i> = t,r,xi,X2,x^. 
We can label various kinds of perturbations according to their transformations under the symmetry group of rotations 
in the 1 — 2 plane. There are three types of decoupled excitations corresponding to spin 2 (scalar channel), spin 1 
(shear channel) and spin (sound channel). 



A quick way to get the first equality is to consider the action of the most general static planar symmetric metrics, vary it, and focus on 
the terms involving second derivatives. Note that with this approach, we have also accounted here for the possible contribution of the 
higher derivative terms in the generalized Gibbons- Hawking term [33| . 
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1. Shear channel 



The shear channel excitations involve hta, hra and h^a with a = 1, 2. Choosing the radial gauge /i^^ = 0, the shear 
channel equations can be reduced to a single equation for Z{r) = qg^^hti + tog^^h^i. At first order in Xgb, Z{r) 
satisfies the equation (below we use the notations introduced in the main text, see p.Sp ) 



Z"(z) + 



4q' 



^2 



q\-z^ + 1) + z4 - 



{z^-\y 



+ 



+ 



Z'{z) 



Z5(g2(^4_l)_^4^)2 



Z(z) 2- 



g2(^ 4^3)_2^4^ 



(Bl) 



Following a similar analysis to that at the beginning of Sec lIII Cl we find that the solution to (jBip which satisfies an 
infalling boundary condition at the horizon z = \ can be written as 



(B2) 



where g is regular at z = 1. In order to find the hydrodynamical poles, it is enough to find g{z) for small values of Co 
and g, which we will assume are of the same order. For this purpose, we introduce a scaled quantity W = and 
expand g{z) as a power series of q. The solution can be readily found to be 



iq 



1 



g(z) = l + 1- 

We thus find near infinity Z{z) can be expanded in 1/z as 

Z{z)k A + Bz-^ + 0{z-^), z^oo, 



where 



= 1 



, . .q , 'HXgbuj 

— — (1 - AXgb) \ . rj. 

Anl uj 4:711 



■ Wq .XcBq ( 1 



AW A W 



2 ^4 



+ o(r) 



i 1 -3A 



GB 



A-kT 



1- A, 



GB 



Carrying out the procedure (|2.6p - (|2.8p one finds that 

Gfl(fc) cx ^ . 

In particular one can show that the poles of Gjiik) solely arise from zeros of A. 

The Dirichlet boundary condition corresponding to ^ = determines the hydrodynamical pole as 
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-iDq^ + 0{q^), 



D = 



Nl_ 
AttT 



(1-3Agb). 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



Note that in the relation (|2.3[) between the diffusion constant D and rj/s, the boundary speed of light c has been set 



to unity (otherwise the right-hand side should be multiplied by c^). Choosing 
so that the boundary speed of light is unity, we find that 

T] 1 



1 — Agb (see equation 



An 



(1-4Agb) + 0(A^b) 



(B9) 



We now need to assume w ~ 0{q^). 
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2. Sound channel 



The sound channel excitations involve hu, hts, hn -\- h22, hm htr, hrs- Choosing the radial gauge /i^^ — 0, the 
sound channel equations can be reduced to a single equation for the variable 



Zs{r) = —g^^hts + 2<?33/i33 - {g^^h22 + ( 1 - 



<7^ drgtt A , o 9^ ^tt 



2 _i "^t. 

LU'^ drgil J ' gii 



(BIO) 



At first order in Xgb, the equation for Zg{z) can be written as (we use the same notation as in the main text) 



=Z':iz) + Z',iz) 



3§^Z^{1 - 5z4) + g2(9 _ ig_24 _^ 15^8) 

z(-l + z4)(_3-|^z4 + ^~2(_i + 3,-4)) 



+ 292^z6(_2 + 3z4) _ g2(_i ^ z4)(-16 + q^z^{-\ + 3z4)) 

z2(-l + z4)2(-3^z4 + ^2(_i + 3^4)) 



+ A, 



4(27|^z8 + 6g2t;;224(_ll + ^4) ^ q4(_;L^ ^ gg^4 _ 27z8)) 



^5(-3f^z4 + 92(_l + 3^4))2 



+ 



^.(^) 



z6(-l + z4)(-3fVz4 + g~2(_l + 3^4))2 



~,6 



~2 ~ 2 

+ g*(g'(7 + z^){z - 3z5)2 ^ 32(4 _ ^iz^ ^ ^5^8)) _ 4l^^4(_^gQ ^ ^32^4 ^ ^2^2(_^o ^ 9^4(3 ^ ^4^^^ 



Again the solution satisfying the infalling boundary condition at the horizon z = \ can be written as 

,{z). 

Defining as above the quantity W — and expanding s(z) in (7, we find that 



1 \ '4]V„ 

^W=(l-;4 



(Bll) 



(B12) 



3W^2^4 _ (1 ^ ^4) -3 + 2z4 + z8 

= 777777^3 —A ^GB- 



iq 



(31^2 _ 2)z4 
T^(z4 - 1) 



z8(3H^2_2) 



z4(3M^2 _ 2) 



XgbW [ 1 



1 \ (-7 + 3(31^2 _ 5)^4) 



4z4(3W^2 _ 2) 



+ 0{q 



(B13) 



The leading asymptotic behavior close to the boundary at infinity is 

Zs{z)^As+B,z-^ + 0{z-^), 

with 



J / 15 
As oc q2(l + Acs) - ^9^^ ( 1 - Y'^^^ 



3tj2 iQXcB 



Again, the hydrodynamical pole is found by setting As = 0, leading to 

t^sound = ±Csq- iTsq'^, 

1 ,r /, , Acs ^ 

c. = ^Arj(i + _), 



2 Ni 



3 47rr 



(1-3Agb). 



(B14) 

(B15) 
(B16) 

(B17) 



By choosing the boundary speed of light to be unity, i.e. Nf^ — a ^ (1 ~ ^f^), we thus find that Cs — and from 



77 1 



(1-4Agb) + 0(A^s) 



s Att 

in agreement with the results obtained from the shear channel and the main text. 



(B18) 
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APPENDIX C: DERIVATION OF (fSTTO]) 



In this appendix we give some details for obtaining gi{z) in equation (j3.19p . Plugging (j3.18p into the equation of 
motion (|3.1ip one finds a fairly complicated ordinary differential equation (ODE) for gi{z). But, by changing variable 
a few times, it reduces to a simpler one. Namely, defining 



1 -4Agb +4Agb^, w = l-it, (Cl) 
we get 

(1 - v){dy{vd,gi + 1)) + 2{vd.,gi + 1) = . (C2) 

Here, we note that — ln(w) is a (singular) solution, as one can also show from more abstract reasoning. In fact, this 
led to our choice of change of variable. Now, we will solve this equation. Defining 

hiiu) ^ {u - l)dugi + I, (C3) 

we have 

uduhi = 2hi, (C4) 

which leads to 

hi = ciu^, (C5) 
where ci is an integration constant. Thus wc find that 

C\U^ — 1 

dugi — ; — = u+1 choosing ci = 1 . (C6) 

u — 1 

Note in order for gi{u) to be nonsingular at the horizon m — 1, we need to choose ci = 1 as we have done above. 
Thus we have 

gi = + u + C2 . (C7) 
We will choose the integration constant C2 so that 51 — > as z oo. This then leads to p.l9p . 



APPENDIX D: STRETCHED HORIZON APPROACH 

In this section, we calculate 77/s for Gauss-Bonnet gravity by extending the stretched horizon approach of @ (see 
also [131). Along the way, we also explicitly show that ry/s is independent of ai and a2 at linear order, as expected 
from the field redefinition argument in Sec lII CI As a spin-off of this work, the framework constructed here allows us 
to consider tipping of the graviton cone in a more abstract way than that presented in Sec lIV Al 



1. Kaluza-Klein reduction 

The stretched horizon calculation of [33] begins with an effective Kaluza-Klein reduction of the AdS black hole 
metric and treating a certain class of off-diagonal metric perturbations as a vector in the reduced geometry. In order 
to develop the effective Maxwell action for h^y, we reduce along the y-direction: 

ds^ = g^.^dx^'dx" + e^P{dy + A^dx^'f . (Dl) 

To construct the theory for a higher curvature action, we need to evaluate the various components of the Riemann 
tensor. This is most efficiently done using an orthonormal frame, i.e. ds^ — tiabE^^E^ , which we can conveniently 
choose as 



= e'^i^i_dx^ with a = i,x,z,f 
Ey = ePidy + A^dxi"), 



(D2) 
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where e"^ are some choice of tetrad components for the reduced metric g^j/, which need not be specified. 
Straightforward calculations then yield the following results: 

Rabcd = Rabcd ~ ■7^^'' {Fa{cFh\d\ ~ FabFcd) 

= [R^\abcd - 2^^'' {Fa[cFb\d] - FabFcd) , 

Rayby = a^ bP - ^ aP^ bP + ^e^" FacFb" 

1 



R, 



abcy 



[RXyby + r'"^''^^'^ ' 
-^eP (^\7cFab + 2S/cpFab + S/bpFac-^apFbc) • (D3) 



Our notation here is such that Rabcd and Vo denote the curvature components and covariant derivative of the four- 
dimensional geometry specified by g^i, . We have also presented the first two curvature components using the notation 
[R^]abcd 'which denotes to the background curvature, i.e. the curvature of the full fivc-dimensional geometry with 
= 0. Hence, for example, [R^]abcy — ^■ 
For later convenience, we also present the components of the Ricci tensor and scalar here: 

Rab — R acb R^ ayb 

= Rab - ^a^bP - ^aP^bP - ^e^" FacFb" 

= [RXb-l^^FacFb^^ 
Ryy = R-'yay - "V^p - {V pf + ^e^^F^ 

I- Jj/j/ 4 ' 

Ray = R\by = "^e" (v^'F^fc + 3 V> F„, 
R = -R"a + R^ y 

= R-2 \/^p - 2 (Vp)2 - ^e^P F"^ 

^ _ le^pp^. (D4) 

2. Curvature-squared theories 

Given the above results, we can begin to apply the stretched horizon approach to the various curvature-squared 
theories considered above. First, we will confirm that for the R^ and R^jyRP" theories rj/s remains unchanged to 
leading order. Hence we begin with the action (|1.2p with — 0: 

1 /• , — /12 



/ - J "^'^ V^ji^ + R + L^ {aiR^ + a2R^.RP'') ) . (D5) 



The background geometry is a planar AdS black hole with metric as in (13.1 



ds' = -f{r)N^df + — + — {dx^+dy^+dz^) , (D6) 



where the event horizon appears at /(r = r+) = 0. (Note that for the present purposes, we do not have to specify 
/(r) in further detail.) We introduce a metric perturbation /i^^ = and perform a Kaluza-Klein (KK) reduction on 
y as above. Then we wish to expand the action (|D5[) to second order in the perturbation. Keeping only the quadratic 
terms, the resulting action is 



^ vcc 



d^x./^e^p{--F^- — 



1 [R^] F' + a, (2 [R^FacFb^ - [R^^ F' - (g-^" V'^ {e'^Fab 



(D7) 
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Now we begin by noting that we are working perturbatively to linear order in ai 2 and that the leading order equation 
of motion for the vector perturbation is: V" {e^''Fab) = 0{ai). As a result, we easily see that the contribution of the 
last term in the above action to the equations of motion will necessarily be 0(a?). Hence this term can be dropped 
in the present analysis. Further, since the background metric (jD6l) will satisfy Einstein's equations to leading order, 
i^^lfiu ~ ~ (4/^^) 5/j,!v + 0(ai). We can make this replacement for the background curvatures appearing in the 0{ai) 
terms in the action, with the result: 



1 



F^\l + 40q!i - 8a2' 



(D8) 



Thus, to linear order, the only effect of these two curvature-squared terms is to change the normalization of the 
effective Maxwell action. The subsequent analysis will be identical to that presented in l6[ with the standard result 
that rj/s = l/47r. 

Next we need to construct the effective action for vector perturbation in Gauss-Bonnet theory (|1.4p . For this 
purpose, we can use the contributions calculated in the above action (|D7p with ai = Xgb/'^ and a2 = ~2Xgb- Next 
we must determine the contribution coming from the Ricmann-squared term. Using the results in (|D3p . we have 



I' 



i^cFab + 2 VcP Fab + Vbp Fac ^^aP Fbc 



(D9) 



where in the end we will substitute 013 = Agb/2. The first term has already been simplified using the cyclic identity, 
R[abc]d = 0. Now the second line above can be simplified by judiciously integrating by parts, applying various identities 
and using the results in (|D3p and (|D4[) . For example, up to total derivatives, we have 



In any event, the final result is 
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The quadratic action for the vector potential arising from the Gauss-Bonnet theory (|1.4p is thus 
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3. Shear viscosity via membrane paradigm 

Next we need to extend the analysis of @ to accommodate the generalized vector action (|D12p which arises in 
Gauss-Bonnet gravity. In particular, we can write the latter in the form 




(D13) 



where the background tensor X necessarily has the following symmetries: 

j^abcd j^[a6][cd] j^cdab (ID14^ 

Though this generalized action will not accommodate the contributions of an arbitrary higher curvature term (e.g., 
compare with (|D7p ). this special form would also apply for generalized Lovelock theories of gravity. Further, in the 
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present case where we are studying a static AdS black brane background, this tensor satisfies the two important 
properties: First, X is diagonal in the index pairs \ah] and [cd], e.g., X*^^^ = 0. Second, all of the components 
of X are nonsingular at the horizon r = r+ when described with frame indices. Alternatively, if the tensor carries 
coordinate indices, the latter can be phrased as saying that all of the components of X^^p" are nonsingular at the 
horizon. 

Given the above framework, one easily extends the analysis of Q. After defining a stretched horizon aX r = th 
(with th > r+ and r/f — r_|_ <C the natural conserved current to consider is 
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(D15) 



where Ua is an outward-pointing radial unit vector. One then simply follows each of the steps appearing in [6] to 
arrive at the following simple result for the effective diffusion constant: 
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For the standard effective Maxwell action with Lagrangian 
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Xt^ 



and then the first expression above reduces to the usual result, first derived in [y]. 

Now we apply this analysis to the specific action (|D12p arising from Gauss-Bonnet gravity. First we must extract 
expressions for the background tensor, which is simplified if we divide up X into three contributions 
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(D17) 



where X'"' and X'-'-^ correspond to the contributions coming from the terms proportional to F'^ and FacFi,'^, respec- 
tively. X<^' captures the remaining contributions. For the action (jD12p . one finds 



yW cd 
^ ab 
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(D18) 



The expression (jD16p for the diffusion constant requires three of the components of X in particular. Using the 
background metric (|D6p and the expressions (jPlSp . one finds that 
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(D19) 



To proceed further, we must explicitly introduce the solution p.2p 
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t2 

2Xgb 



1 - A 1 - 4A 
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(D20) 



for the black brane in the Gauss-Bonnet theory. Recall that the temperature (|3.4p is given by T = iVj r+/ ttL^. Further 
implementing the KK reduction (|Dip on this background (|D6|) yields 
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(D21) 
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Given these results, the prefactor in (jD16p reduces to 



r— r-|- 



(D22) 



We note that the second factor in Xtr has a particularly simple form: 1 — 2Ags-t-/ — \ — 4Agb I 1 i 



Then the integral in (|D16p is evaluated as 

i~)dr 
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GB 



Combining the results in (|D22p and (|D23p then yields 



4r+ 2Ag_b 47r 1 



(D23) 



(D24) 



where c — N^/a is the boundary speed of light, with a defined in p.3p . Hence we recover the expected result for the 
Gauss-Bonnet theory: 
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4. Graviton cone revisited 



Given that the background tensor X is expressed in terms of curvatures of the background spacetime (see (jD17p 
and (jPlSp ). we should be able to express the effective "null" cone of the gravitons in terms of these curvatures. 

In Sec lIV Al we considered the scalar channel which corresponds to a perturbation hx^ with dependence on t, r and 
z. In the present notation, this is precisely an excitation of the vector component A^- In a high-frequency or WKB 
limit, we write 

Ax = e'' " 0en , (D26) 

where the first factor is the rapidly varying phase and 0eri is the slowly modulated envelope function. The coordinate 
dependence of the scalar channel also requires that kx — — ky. For these modes, the efi'ective action (jD13p reduces 
to 




(D27) 



From this action, we can readily derive the full equations of motion, however, we do not need these here. In the 
high-frequency limit, the graviton cone is given by 

= X"^''^ ka kb . (D28) 

Now as indicated above, we use (jPTTl and (IdTsI) to express this result in terms of the background curvatures. Hence 
the effective metric defining the graviton cone can be written as 

~^2XgbL^ [[RXy'^ + [RXx"' - 2 [R%x^' Sa") , (D29) 

where we have implicitly assumed that a and b only take values in {i, f , z}. We have also canceled certain terms using 
Rx^ = Ry^ for the backgrounds of interest here. In the first line, the correction term is proportional to the Einstein 
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tensor, a result that is reminiscent of that in [52, l53|. Their results for the characteristic hypersurfaces of Gauss- 
Bonnet gravity do not include the nontrivial contribution in the second line above. We do not entirely understand the 
source of this discrepancy but note that the analysis of [sl, [s^ uses complementary techniques to ours. At least in the 
context of small \gb, the additional terms in (|D29p have an important consequence. That is, using the equations of 
motion for the background geometry, the results of 52] would have predicted that the deviation of the graviton cone 
from the standard light cone only occurs at 0{\^g). However, our results in Sec lIVI indicate that there is nontrivial 
result at 0{\gb)- The additional terms appearing in the second line of (ID29|) must be responsible for this effect. 
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